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Abstract 

We consider the inverse problem of determining the time dependent magnetic 
field of the Schrodinger equation in a bounded open subset of M n , n > 1, from a 
<D \ finite number of Neumann data, when the boundary measurement is taken on an 
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appropriate open subset of the boundary. We prove the Lispchitz stability of the 
magnetic potential in the Coulomb gauge class by n times changing initial value 
suitably. 
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1 Introduction 

1.1 Statement of the problem 

Let T > 0, n G N*, and let Q be a bounded open subset of M n with C 2 -boundary 
T. We consider the time-dependent Hamiltonian H a (t) : = (iV + x(^) a ) 2 associated to 
the nondivergent magnetic vector potential x(i)a(x), where x is a smooth real valued 
function on [0, T] and a G H 1 (f2) n is bounded, together with the related Schrodinger 
equation, 

-iu'(t, x) + H a (t)u(t, x) = 0, (t, x) G Qfe := (0, T) x ft 
u(t,x)=0, (t, x) G £+ := (0, T] )x r (1.1) 

u(0, x) = uq(x), x G Q, 

for some suitable data u . Here and throughout all this text, u'(t, x) stands for d t u(t, x) 
and Vu(t, x) := (d Xl u(t, x), . . . , d Xn u(t, x)) is the gradient of u(t, .) at x — (xi, X2, ■ ■ ■ , a; n ) G 
Q. As follows from Remark 12. II below, we may assume in the sequel with no loss of gen- 
erality that the magnetic potential vector a is real-valued. 
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Let v denote the unit outward normal vector to T and set d v v := Vv.u. Then, r + 
denoting an open subset of T satisfying an appropriate geometrical condition we shall 
make precise further, we aim to retrieve a = (aj(x))i<j< n , x G f2, in (11. ip . from the 
extra data d u d^u\ / 0T ) xr+! k — 1, 2, by n times changing initial value u suitably. Hence 
we investigate the problem to know whether a finite number of partial Neumann data 
of fll.ip (in absence of any Cauchy lateral data information given by the Dirichlet to 
Neumann map, denoted by "DN map" for short in the sequel, A a , associated to a) 
determines uniquely the magnetic potential a in the Coulomb gauge class (i.e. the class 
of divergence free vectors in Q). 

1.2 Existing papers 

Numerous papers establishing the uniqueness of inverse problems coefficients from the 
DN map (or scattering information) have actually been published over the last years. 
In the particular case of the magnetic Schrodinger equation, it is noted in [1] that the 
DN map is invariant under the gauge transformation of the magnetic potential, i.e. 
A a +v^ — A a when ip G C 1 (f2) is such that fy? = 0. Therefore the magnetic potential 
cannot be uniquely determined from the DN map (we can at best expect uniqueness 
modulo a gauge transform of a from A a ). However the magnetic field da, where da is 
the exterior derivative of a interpreted as the 1-form X^j=i a jd£j, is preserved. If n = 3 
then da corresponds to curl a. Conversely it is shown in [2] for multiconnected domains, 
that if the DN maps A a and A a are gauge equivalent (i.e. e~ 1 ^'A a e 1 ^ = A a ) then a and 
a are gauge equivalent too. 

Actually Z. Sun proved in [3] that the DN map determines the magnetic field provided 
a is small in an appropriate class. In jl] the smallness assumption was removed for 
C°° magnetic potentials. The regularity assumption on a was weakened to C 1 in [5] 
and to Dini continuous in [6], and the uniqueness result was extended to some less 
regular but small potentials in [7]. Recently in j8], M. Bellassoued and M. Choulli 
proved that the magnetic field depends stably on the dynamical DN map. In [9J, G. 
Eskin considered the inverse boundary value problem for the Schrodinger equations with 
electromagnetic potentials, in domains with several obstacles. He proved the uniqueness 
modulo a gauge transform of the recovery of the potentials from the DN map, under 
geometrical conditions on the obstacles. In [10], M. Salo reconstructed the magnetic field 
from the DN map using semiclassical pseudodifferential calculus and the construction 
of complex geometrical optics solutions. All the above cited papers considered time 
independent magnetic potentials. The uniqueness in the determination from the DN 
map, of time dependent magnetic potentials appearing in a Schrodinger equation (in a 
domain with obstacles), was proved by G. Eskin in [IT]. The main ingredient in his proof 
is the construction of geometric optics solutions. As far as we know, this is the only 
existing paper dealing with the determination of a time dependent magnetic potential 
in a Schrodinger equation. 

All the above mentioned results were obtained with the full data, i.e. measurements 
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of the DN map are made on the whole boundary. The uniqueness problem by a local DN 
map was solved by D. Dos Santos Ferreira, C. E. Kenig, J. Sjostrand and G. Uhlmann 
in [12]. Similarly it was shown in [13] that the magnetic field depends stably on the DN 
map measured on any subboundary T which is slightly larger than half the boundary. 
This result was extended in [J3] to arbitrary small T provided the magnetic potential 
is known near the boundary. 

Notice that infinitely many observations were required in all the above cited results. 
To our knowledge (and despite of the fact that this seems more suited to numerical 
applications) there is no result available in the mathematical literature on the recovery 
of a magnetic potential appearing in a Schrodinger equation, from a finite number of 
boundary measurements. Nevertheless this is not the case for external electric potentials. 
Indeed, the problem of stability in determining the time independent electric potential in 
a Schrodinger equation from a single boundary measurement was treated by L. Baudouin 
and J. -P. Puel in [15]. This result was improved by A. Mercado, A. Osses and L. Rosier 
in [TH]. In the two above mentioned papers, the main assumption is that the part of 
the boundary where the measurement is made satisfies a geometric condition related to 
geometric optics condition insuring observability. This geometric condition was relaxed 
in [T7J under the assumption that the potential is known near the boundary. 

In the present article we prove Lipschitz stability in the recovery of the time depen- 
dent magnetic potential appearing in the Schrodinger equation, from a finite number 
of observations d v d^u, k = 1,2, measured on a subboundary for different choices of uq 
in fll.ip . by a method based essentially on an appropriate Carleman estimate. We refer 
to [18] , [15] and [19] for actual examples of this type of inequalities for the Schrodinger 
equation. The original idea of using a Carleman estimate to solve inverse problems goes 
back to the pioneering paper [20] by A. L. Bugkheim and M. V. Klibanov. This tech- 
nique has then been widely and successfully used by numerous authors (see e.g. [2Tj . 
[22], [23], [23], [15], [25], [26], [27], [28], [29], [30], [3]], [32], [33] and references therein) in 
the study of inverse wave propagation, elasticity, or parabolic problems. However, due 
to the presence of time dependent coefficients (involving the magnetic potential vector a 
we aim to retrieve) for zero and first order space derivatives in the expression of H a (t), 
the solution to the inverse problem we address in this text cannot be directly adapted 
from the above references. 

1.3 Main results 

In this section we state the main result of this article and briefly comment on it. 

Choose a £ H 1 (0) n nH div o(0;l), where H div0 (fi;M) := {a £ L°°(fi;M n ), V.a = 0} is 
the space of real valued and bounded magnetic potentials vectors in the Coulomb gauge 
class, and define the set of "admissible potential vectors" as 

A(a , M) := {a £ H 1 ^)" n H div0 (fi; R), ||a|| Loo{n) n < M and a(<r) = a ((x) a.e. o £ T}. 
By selecting a in A(ao, M) we enjoin fixed value to a on the boundary, which is nothing 
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else but the measurement on T of the magnetic potential we want to determine. Similar 
(or even stronger) "compatibility conditions" imposed on inverse problems coefficients 
have already been used in various contexts in e.g. [31], [32], [11] or [T7j . 

The main result of this paper is the following global stability estimate for magnetic 
potential vectors in A(a , M). 

Theorem 1.1. Let T > 0, n, Q and T be the same as in §1.11 Let T + C T ful- 
fill the geometrical condition of Assumption [X71 let a G H 1 (f2) n fl Hdi V o(^;^); an d 
let x G C 3 ([0, T]; M) be such that x(0) = and x'(0) ^ 0. Pick n functions Uqj G 
j^max(6,n/2+i+e)^_ M.) ; j — 1, . . . , n, for some e > 0, satisfying: 

det DU (x) ^ 0, x = (xi, . . . ,x n ) G Q, where DU {x) := {d X] u 0)i (x))i<ij< n . (1.2) 

Let a (resp. a.) be in A(a , M) , and let Uj (resp. Uj) denote the C°([0, T]; Hj^f^nH 2 ^) )- 
solution to ( 11.11) (Vesp. ( II. II) where H a (t) is replaced by H a (t) ) with initial condition Uoj, 
j = 1, . . . , n. Then there exists a constant C > 0, depending only on T, f2, T + , M, x, and 
{wo,j}j=i> sitc/i i/iai we /iai>e: 

n 

||a - a|| L 2 (n) „ < {\\ d Muj - Mj)||L2( ,T;r+) + R^tOj - Ui)||tf( 0l T;r+)) ■ 

i=i 

This immediately entails the: 

Corollary 1.2. Under the conditions of Theorem the following implication holds 
true for every a and a in A(a , M): 

(uj(t, x) = Uj(t, x), j — 1, . . . , n, a.e. (t, x) G (0, T) x T + ) (a = a) . 

These results suggest several comments. 

1. We emphasize that the stability result of Theorem II . 1 1 only requires a finite number 
(directly proportional to the dimension of Q) of observations on an appropriate 
subboundary, for a finite span-time. This result actually involves less measure- 
ments than the other reconstruction methods for magnetic potential vectors (or 
magnetic fields) known so far. 

2. Although the inverse problem examined in [T5], of determining the external po- 
tential q appearing in the Schrodinger equation — \u' + Am + qu = in Qt, where 
Am := Y^j=i ®x u 1S the Laplacian of u, from the single measurement might 
seem very similar to the one addressed in this paper, there is an additional major 
mathematical difficulty when dealing with the magnetic potential. This comes 
from the presence of first order spatial differential terms of the form a.Vu in the 
effective Hamiltonian H a {t). 
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3. If n = 3, the equations in (11. ip model (in a "natural" system of units where 
the various physical constants are taken equal to 1) the evolution of the wave 
function u of a charged particle subject to the action of the magnetic field b(t, x) : = 
x(t) curl a(x), starting from the state u . Notice that in classical physics, only 
b has a physical meaning, and a is a mathematical tool only. Moreover, in light 
of [3"5] [Remark 2.1], one can always assume (at least for a sufficiently smooth T) 
that a satisfies the Coulomb gauge condition V.a = 0. As already mentioned in 
§1.11 there is an obstruction to the uniqueness of a, since the magnetic potential 
appearing in the Schrodinger equation can at best be recovered modulo a gauge 
transform. The Coulomb gauge condition imposed on a actually eliminates this 
indetermination. 

4. In the particular case where — sin(cjt) for some to > 0, (jl.ip describes the 
evolution of the wave function u of a (periodically) laser pulsed charged particle 
(see e.g. [36][Chap XVIII,§1.4.2.2]). 

5. The Lipschitz stability inequality stated in Theorem 11.11 is not only interesting 
from the mathematical point of view, but may also be a very useful tool in view 
of numerical simulations (see [37] where this point is discussed). 

1.4 Contents 

The paper is organized as follows. In section [2] we study the solutions to non autonomous 
magnetic Schrodinger equations. Existence, uniqueness and regularity results are stated 
in §2.41 and we establish in §2.4.31 that the charge and the energy of these systems 
remain uniformly bounded in the course of time. The main properties of the associated 
magnetic Hamiltonians, needed in the proofs, are collected in §2.21 and §2.31 Section [3] 
is devoted to establishing the stability inequality stated in Theorem 11.11 The inverse 
problem of determining the magnetic potential vector from partial lateral Neumann data 
is discussed in §3.11 The strategy used essentially relies on a global Carleman estimate 
for magnetic Schrodinger equations, given in §3.21 The proof of the global stability 
inequality is given in §3.31 

2 Non autonomous magnetic Schrodinger equations 

The system under study is modeled by non autonomous magnetic Schrodinger equations 
which do not fit the usual requirements (see e.g. [38] [Chap. 5, Assumption (12.1)] 
imposing time independent coefficients to (non zero) order space derivatives) of the 
classical existence, uniqueness or regularity results in L 2 (Qj), found in the mathematical 
P.D.E. literature. Furthermore, the strategy used in Section [3] requires that the system 
(II. ip be differentiated twice w.r.t. t, and, roughly speaking, then solved in a wider space 
of the form L 2 (0, T; wh ere W denotes the dual space of Hj(fi) or Hj(fi) n H 2 (ft). 



5 



This statement is made precise in Lemma T2.41 The corresponding existence, uniqueness 
and regularity results on the solutions to these problems are collected in Lemma 12. 2\ 
while the question of their equivalence in the above various functional spaces is treated 
by Lemma 12.31 Finally, some a priori estimates used in the derivation of the stability 
equality of Theorem 11.11 are derived in Lemma 12.51 and Proposition 12.61 

2.1 Notations 

In this subsection we introduce some basic notations used throughout the article. Let 
Xi, X 2 be two separable Hilbert spaces. We denote by B(X 1 ,X 2 ) the class of linear 
bounded operators T : X\ — > X 2 . Let (a, b) be an open set of R. If the injection 
X\ X 2 is continuous and X\ is dense in X 2 , we define W(a, b; Xi, X 2 ) := {u; u G 
L 2 (a, b;Xi), u' G L 2 (a, b;X 2 )}, which, endowed with the norm \\u\\w '■= ( 1 1 ^ 1 1 l 2 (a b-Xx) + 
ll M 1lL 2 (a 6-x 2 )) 1 ^ 2 ' * s a Hilbert space. More generally we put W m (a, b; Xi, X 2 ) := {u; u G 
H m (a, b;Xi), u' G H m (a, b; X 2 )} for every m G N*, where H m (a, b; X\) denotes the usual 
m th -order Sobolev space of Xrvalued functions. For the sake of convenience we set 
HV^Xi) := L 2 (a,6;X!) in such a way that W°(a, b] X x , X 2 ) := W(a,b;X 1 ,X 2 ). 
If X x = X 2 = X we write B{X) instead of B(X,X), and W m (a,b;X) instead of 
W m (a,b;X,X), m G N. 



2.2 Magnetic Schrodinger operators 

Let a G C°([0, T]; Hdivo(^l We consider the linear self-adjoint operator H(t), 

t G [0,T] in "H : = L 2 (fi), associated with the closed, densely defined and positive 
sesquilinear form 

h(t; u, v) := ((iV + a(t, x))u, (iV + a(f , x))v) , (2.1) 

with a domain independent of t, dom h(t) = dom h(0) = Hj(fi). Here (., .)o denotes the 
standard scalar product in Hq. The space dom h(0) endowed with the scalar product 
(u,v)x := ((-A + ly^u, (—A + l) l/2 v) (or equivalently (u,v)i = (u,v) + (Vu,Vv) ) 
is a Hilbert space denoted by Tii, and we have 

(H(t)u, v)o = h(t; u,v), u G dom H(t), v G Hi. (2.2) 

The boundary T being C 2 , we actually know from [39] [Chap. 2] that 

dom H{t) = dom (-A) = Hj(fi) n H 2 (fi). (2.3) 

We call the dual space of Hi, that is to say the vector space of continuous conjugate 
linear forms on Hi- For any u G Ho, the functional v i— >■ (u,v)q belongs to since 
|(^,f)o| < ll n llolK'||o < IMIolMli) an d we can also regard Hq as a subspace of H-\. 
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Hence "Hi C "Ho C where the symbol C means a topological embedding, and H(t), 
t G [0, T], can be extended into an operator mapping "Hi into H-i since 

\h(t;u,v)\ < Ct\\u\\i\\v\\i, u,v G Hi, 

where 

C T := (1 + nAl) 1/2 , A := ||a||co([ 0) T];L°°(n)«) = sup \\a(t, .)||L°°(n)*- (2.4) 

*e[o,T] 

Let us denote by (■, ^-i^ the dual pairing between H-\ and H\. This pairing is linear 
in the first and conjugate linear in the second argument. In other words, the embedding 
"Ho C H-i means that (g, "0)-i,i — (flS^o f° r ah 9 ^ Ho an d ^ £ an d the mapping 
-fT(t) : "Hi — > H-i is defined so that (H(t)u, v )-i t i = h(t;u,v) for all u,v G "Hi. 

The space dom (—A), endowed with the scalar product (u, v) 2 ■= ((—A + l)u, (—A + 
l)u) , is an Hilbert space denoted by H 2 . For all u G H 2 , we have 

if (t)u = (iV + a(t)) 2 w = (-A + 2iRe (a(t)) .V + |a(t) | 2 )w, (2.5) 

from ( I2.2p -( l2~3"j) and the Coulomb gauge condition V.a = 0, which entails 

\\H(t)u\\ < \\Au\\ + 2A \\ Vtt|| + ri^olkllo < 2C|||m|| 2 . (2.6) 

We call H- 2 the dual space of H 2 . Similarly, we have H2 C C Ho C H-i C H_2, and 
we deduce from (12.61) that H (t) can be extended into an operator mapping Ho into H~ 2 . 
Denoting (., .)_2,2 the dual pairing between H- 2 and H 2 , the mapping H(t) : Ho — >■ %_ 2 
is defined by 

(H{t)u,v)_ 2>2 := (u,H{t)v} , u G H , u G ft 2 . 

Remark 2.1. In light of (12. 5 p we may assume in the sequel without limiting the gen- 
erality of the foregoing that the magnetic vector potential a is real- valued. 



2.3 Time-dependent magnetic Hamiltonians 

In this section we examine the dependence of H{t) w.r.t. t and establish some of the 
properties of its derivatives which will play a crucial role in this paper. 

For j = 1, 2, 3 we assume that a G C J ([0, T]; H div0 (f2; K)) and introduce the operator 

B 3 (t) : = a W(f).(iV + a(t)), a^(t) := 

with domain Hi, acting and densely defined in Ho- In light of the divergence free 
condition imposed on a, we have 

(BjWp, q) = (p, B 3 (t)q) Q , p, q G Hi, (2.7) 
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thus Bj(t) is symmetric in "Ho- Further, for all u G Hi, we have 

WBjfyuWo < ^||(iV + a(t))«|| < ^.j-(7t||u||i, (2.8) 
where Ct is the same as in (12.41) . and 

■■= ||a||cj([o.r];L°°(n)") = sup ||a (j) (t)|| LO o (n) n, j = 1, 2, 3. (2.9) 

te[o,T] 

Hence Bj(t), j = 1, 2, 3, can be extended from % Q into H-i, by setting: 

u)_ ltl := (u, 5,^)^)0, u G H , u G Hi. (2.10) 

For all £ G [0, T], we then define the three following operators 

H w {t) := 2Bi(t), # (2) (t) := 2[5 2 (t) + a'(t) 2 ] and H {3) (t) := 2[B 3 (t) + 3a'(t).a"(t)] 

(2.11) 

with domain Tii- In light of (12. 7p the operators H^'{t), j = 1, 2, 3, are symmetric in "Ho 
for each i G [0, T], and, due to (12. 8p . there exists moreover a constant lj > 0, depending 
only on Ct and {^4fc}^. =0 , such that we have: 

||# (i) (0 u llo <^'N|i, ueHi. (2.12) 

Hence H^\t) can be extended from "Ho into H_i as (H^\t)u,v)^i t \ = (u, H^\t)v )o 
whenever w G "Ho and v G "Hi- Moreover for a in C J ([0, T]; H div0 (f2; R)), j = 1,2,3, 
the mapping t H- u, v), for u,v eHi fixed, is j times continuously different iable on 
[0,T], and it holds true that 

^-h(t;u,v) = h®{t;u,v) = (H®(t)u,v) = (u,H^(t)v) , J = 1,2,3. (2.13) 

When convenient, we shall write in the sequel H'{t) for H^(t), H"(t) instead of H^(t), 
and H"\t) for H^\t). 

2.4 Solution to non-autonomous Schrodinger equations 

Throughout this section we assume that a G C 1 ([0, T}; Hdi v o(^; R))- 

2.4.1 Solutions in L 2 (0, T; H-j), j = 0, 1, 2. 

Let / be in L 2 (0, T; %_./) for j = 0, 1, 2. A solution to the Schrodinger equation 

- vj)' + F(t)^ = / in L 2 (0, T; H-j), (2.14) 
is a function ip G L 2 (0, T; %2-j) satisfying for every v EHj: 

-i±(m,v)o+{H{tMt),v)-jj = (/(<)> «>-iJ in <T(0,T)', (2.15) 
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where (., stands for (., .)o in the particular case where j = 0. 

Let us now introduce the following notations that will be used in the remaining of 
this section. For all v G C§°(fi) = Cg°(fi; C) and tp G C°°(0, T) = C%°(0, T; C) we set 

(v <8 ip)(t, x) := v(x)tp(t), t G (0, T), 160, 

and define the space Cg°(fi) <8 Cg°(0,r) := {v <8 </?, t> G CS°(fi), ^ G Cg°(0,T)}. We 
recall that 

Cg°(f2) (8) C£°(0,T) is dense in C£°(Q+) = C~(Q+;C), Q+ := (0,T) x Q. (2.16) 

Let V be solution to f[2TT3j) . For all u G C£°(fi) and (/? G C °°(0, T; E), it follows from 
(EUSD that 

i/ (iP(t),v) <p'(t)dt= [ (g(t),v)-jj<p(t)dt, (2.17) 
Jo Jo 

where #(t) := /(i) - H(t)%j)(t) for a.e. i G (0,T). Evidently 3 G L 2 (0, T;%_j), and 
fETTj) yields 

i/ ^{t),{v®^)'{t)) dt = [ (g{t),{v® V )(t))- jtj dt, (2.18) 
Jo Jo 

for all u (8) <p G Cg°(fi) ® Cg°(0,T). Since if) G L 2 (0,T;% ), £ e L 2 (0, T; "H-j), and 
Cg°(fi) ® Cg°(0,T) C Cg°(Q+), (CTD can be rewritten as 

i(if}, (v (8) v) , )cg»(g+)' ) cg=(o+) = (9(t),v® <p)cf(q+)',c%>(q+)> c o°( fi ) ® c o°(°> T )- 

From this and (I2.16P then follows that i^' = g in C x3 ((5y) / . Taking into account 
that G L 2 (0,T;^_ y ) we obtain that V' e L 2 (0, T; and thus that fl27H|) holds 
true. Henceforth any solution to (I2.14p belongs to W(0,T;'H2-j,'H-j) and (I2.15P can 
be rewritten as 

-iW(t),v)- jd + (H(t)m,v)- hJ = (f(t),v)^j in C~(0,T)', u G ^. (2.19) 

As a consequence a solution to the Schrodinger equation (12. 14ft is a function ip G 
W(0,T;H 2 - j ,'H- j ) satisfying f[2TT9l) . 

2.4.2 Existence and uniqueness results 

For j = 0, 1, 2, we consider the Cauchy problem 

-h// + H(t)4> = / in L 2 (0, T; 2 
where ^o G and / G L 2 (0, T; 
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Notice that the second line in (I2.20p makes sense since ip(0) is well defined in Hi-j for 
every j = 0, 1, 2. Indeed we have W(0, T; H 2 -j, H-j) ^ C°([0, T]; Hi-j), j = 0, 1, 2, as 
can be seen from [36][XVIH.§l.(1.61)(iii)] for j = 0, and from [36] [X VI 1 1. §1. Theorem 1] 
for j = 1, the case j = 2 being a direct consequence of the imbedding W(0, T; "H 2 , "H ) <^-> 
C°([0, T}; Hi) and the fact that (1 — A) -1 maps W(0, T; Ho, H-2) one-to-one onto 

w(o,T;H 2 ,h ). 

Further, with reference to [TTj [Lemma 2.1] for j = 0, and to (38] [Chap. 3, Theorem 
10.1 & Remark 10.2] (see also [36] [XVIII. §7, Theorem 1 & Remark 3]) for j = 1, we 
recall the following existence and uniqueness result: 

Lemma 2.2. Let T, n and Q be as in let a £ C 1 ([0, T]; Hdi v o(^; one? /ix 

j = 0,1. T/ien /or a// ^ £ H 2 _j and f £ W(0, T; Ho, H-j), there exists a unique 
solution ip £ C°([0,T};H 2 - j) C) T}; H-j) to (12^0]) . 

In light of Lemma 12.21 it is thus equivalent to solve (I2.20p with either j = or j = 1 
when ip £ H 2 and / £ W(0, T; Ho)- Moreover, due to the following uniqueness result, 
the same is true for either j = 1 or j = 2 provided ip £ Hi and / £ W(0, T; Ho, H-\): 

Lemma 2.3. Let T, n and Q be the same as in §i. 1\ Then, for all f £ L 2 (0, T;H- 2 ) 
and ipo £ H-i, there exists at most one solution to the problem ( 12 . 20 j) for j = 2. 

Proof. It is enough to prove that the only solution %p to the problem (I2.20p with / = 
and ipo = is identically zero. Actually for all r £ (0,T), the first line of (I2.20p with 
/ = states that we have (i//(t), i>)-2,2 — — i(0( r )> H(r)v) f° r an y u G In the 
particular case where v = R(t)iP(t), with -R(t) := (1 + if(r)) -1 , this yields 

(^)^(rMr)>- ai2 = i(||fl(r) 1 ^(r)||g-||V'(r)||g). (2.21) 

In light of (I2.13p . r i-> R(r)ip is differentiable in (0,T) for each £ "H , and we get 
4-R(t)<p = R(t)H'(t)R(t)(p by standard computations. Therefore r H- (R(t)iJ}(t), <p) = 
(^}(r), R(r)(p)o is differentiable in (0,T) as well, and we find that 

^( J R(r)^(r),^) = (^(t),R(t)v)-2,2 + (4>(t),R(t)H'(t)R(t)v)o 

= (^(r), i?(r)^)_ 2i2 + {R{t)H'{t)R{t)^{t), <p) Q . (2.22) 

Recalling that i?(r) is extended into an operator mapping Hq into H-2, defined by 
(H(t)u, u)-2,2 = ( M ) H{t)v)q for all u EHo and i> £ "H 2 , the resolvent -R(r) can therefore 
be extended into an operator mapping H-2 into "Ho, according to the identity 

(R(t)u,v) := (u, R(t)v) -2,2, u £ H-2, v £ Ho- 

This, together with ( 12.22ft proves that (0, T) 9 r H- R{r)ip{r) is differentiable in 
"Ho, with ±R{t)iI){t) = R{t)^'{t) + R{t)H'(t)R(t)iIj{t). As a consequence, r ^ 
|| J R(r) 1 / 2 ^(r)||2 = (^(t),R(t)iP(t))o is differentiable on (0,T), and we have 

±\\R( T y/ 2 i;(T) \\l = (^(r), R{t)H>{t)R{t)4,{t))o + 2Re «^( T ), i?(r)^(r))_ 2 , 2 ) . 
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The second term in the r.h.s. of the above identity being zero according to (I2.2ip . we 
end up getting that 

^WRiT^WWl = (R(t)iP(t),H'(t)R(t)i/)(t)) . (2.23) 

Bearing in mind that H'{r) = 2B 1 (r) and using (12.91) we get that 

WH'^R^rf^W < 2A U r G (0,T), (2.24) 

since WB^Riry^uWo = ||a'(t).(iV + a^))^) 1 / 2 ^^ < A || (iV + a^))^) 1 / 2 ^^ for 
all u G Hq, and 



(iV + ^Rirf^uWl = WHirY^RirY^uWl < 



u 



Putting and ([221 together, we obtain that ^||i?(r) 1 / 2 ^(r) || 2 < 2Ai \\R{t) 1 / 2 ^ 

for every r G (0, T), hence 



\\R{t) 1/2 m\\i<2Ai / \\R(Ty/^(T)UdT, te[0,T], 

Jo 

by integrating w.r.t. r over [0, t] and using the fact that ^(0) = 0. Therefore i?(t) 1/2 ^(t) = 
for all t G (0, T) by Gronwall inequality, which yields ip = and proves the result. □ 

In light of section §2.31 it is possible to differentiate (at least in a formal way) the 
Cauchy problem (12.201) w.r.t. t. This is made precise with the following: 

Lemma 2.4. Let T, n, Q and a be as in Lemma \2.2i, let f G 1^(0, T; Hq), i/)q G 
H 2 , and let $ denote the C°([0,T]; ft 2 ) n C 1 ( [0, T] ; T-L )-solution to fl2^0|) . Then \\)" G 
L 2 (0, T; %~2), and ip' is solution to the system: 

-if + H(t)if)' = f - H'(t)i/) in L 2 (0, T; U- 2 ) (t) 

^'(o) = -i(#(oM)-/(o)). 1 ^ 

Proof. The existence and uniqueness of if) G C°([0, T]; H 2 ) H C 1 ([0, T]; %q) being guar- 
anteed by Lemma |2T2"1 j = 0, the mapping t \-> (H(t)ip(t),v) = (if>(t), H (t)v) is thus 
continuously differentiable in [0,T] for every v G and it holds true that: 

^ t (H(t)if)(t), v) = (if)'(t), H(t)v) + (i[>(t), ff'(t)v) = + U>-2,2. 

Further we have t >->■ (f,v) G W(0,T;C) by assumption, hence the same is true for 
t ^ (ip'(t),v) from (EZOD , and the result follows. □ 
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2.4.3 Uniform time boundedness of charge and energy revisited 

We establish in Lemma 1231 1 hat the charge and the energy of the system driven by (I2.20p 
remain uniformly bounded in the course of time over [0, T]. These estimates, which are 
reminiscent of [3"o] [Chap. XVIII, Formula (7.17)] or [To] [Lemma 7], are obtained here in 
the more general L 2 (0, T; "H.^-framework for (I2.20p . They are further extended to the 
cases of the first and second time derivatives of the above quantities in Proposition 12.61 
This result is one of the main ingredients in the derivation of the Carleman estimate of 
Proposition 13.31 

Lemma 2.5. Let T, n, Q and a be as in Lemma \2.2 . Then, there is a constant 
Co > 0, depending only on T, Ao, A\ and Q, such that for all ipo G Hi and all 
f G W(0,T;^ ,H_i), the C°([0,T];fti) n C^QO, T]; U^) -solution rj) to fl2^U|) . satis- 
fies simultaneously: 

\\m\\o < e T/2 (U \\ + II/IIlw^o)), t G [0,T], (2.26) 

and 

\\m\\i < co(ll^o||i + \\f\\w(o,T;H ,n^)), t G [0,T]. (2.27) 

Proof. We have -i(^ (t),^(t )) + h{t-^{t)^{t)) = (f(t),^(t)} for all t G [0,T], by 
multiplying (I2.20p by ip(t,x) and integrating w.r.t. x over Q. Selecting the imaginary 
part in the obtained expression then yields Re ({if)'(t), i/j(t))o) = —lm((f(t),i/}(t))o), or 
equivalents ^||V>(*)||o = - 2Im ((/(*)> ^(*))o). This entails 

\\m\\l<\\Ml+\\f\\h { o,T-,H )+ f \M*)\\ods, te[o,T], 

Jo 

and (I2.26P follows immediately from this and Gronwall inequality. 

We turn now to proving fl2T27p . The set W(0, T; Ho) being dense in W(0, T; H , H-i), 
we consider a sequence (f n ) n of W(0, T; Ho) satisfying ||/— fn\\w(o,T;H ,H-i) — ^ as goes 
to infinity. Similarly we pick a sequence (ip n ,o)n of Hi such that lim^oo ||^ n ,o — V'olli = 0, 
and, for each n G N, we call the C°([0, T]; "H 2 )nC 1 ([0, T]; "H )-solution to (I2T201) . where 
(/,^o) is replaced by (f n ,^n,o)- 

f -W n + H(t)i> n = f n in L 2 (0, T; H ) 

\ 1p n {0) = -0n,O- 

From f l2T2Tjp . ([2726]) and fl2728|) then follows for every n G N that 

||^(t) - ^(t)||o < e T / 2 (||^0 - V'n.oHo + ||/ " fnh*(0,T;H )), t G [0, T] 

Further, multiplying (12.281) by ip' n (t,x), integrating w.r.t. x over Q, and selecting the 
real part of the deduced result, entails Hje((H(t)i(} n (t),i(}' n (t))o) = Re ((f n (t), ip' n (t))o) 
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(2.28) 



(2.29) 



for all t G [0,T]. Hence ±h{t; ^ n (t), if> n {t)) = h'(t; ^ n (t), ^ n (t)) + 2Re {{f n (t), #,(t))o), 
which yields 

/i(t; if) n (t),i/) n (t)) - h(0; -0n,o, V'n.o) 
= | / l '( S ;^( S ),^( S ))d S + 2Re(|| (/^^^lo) ds, te[0,T]. (2.30) 

Moreover, as / *(/„(s), <(s)) ds = (/„(*), ^ n (*))o - (/n(0), ^ n ,o)o - / *(/£(s), Vn(s))ods, 
we have 



< ||/„(t)||-l|WOIIl + ll/»(0)ll-l|l^n.o|| 



+ / ||£(*)||-i|hM*)||id* 

Jo 

< Ci(l + e~ 1 )\\fn\\w(!a,T;'Ho,H-i) + II^MIl 

+e||^(*)lli+ / ll^n(s)||?d S , (2.31) 
Jo 

for every e G (0, 1), the constant C\ > depending only on f2. Here we used the fact 
that W^T-U^U-!) C°([0,r|;?{_i) (see [36] [Chap. XVIII §1, Formula (1.61) (iii)]) . 
Further, putting (12.261) and (I2.30[) - (I2.31[) together, and taking into account that h is %i- 
coercive w.r.t. Ho, we may choose e > small enough so that we have 

Un(t)\\l <°2 (|hMI?+ ||/n|| jf ||^(s)||?dsj , te [0,T], 

for some constant c 2 > depending only on T, ^4o, .Ai and fl. Therefore we end up 
getting 

||^„(*)||l < Codl^oll? + ||/n||^(0,r;« 0| «- 1 ))' 1 E M' ^ 

from the Gronwall inequality, where cq fulfills the conditions of (12.271) . 

From (12.281) and ( 12. 32ft then follows that {ip n )n is a bounded sequence in L°°(0, T; Hi), 
and ("0n)n is bounded in L°°(0, T; "H_i). Moreover, for every £ G [0,T], (VVi(£))n con- 
verges strongly to ip{t) in "H by (12.291) . hence weakly in %\ according to [39J [Proposition 
1.3. 14(i)] . Therefore we have H^WIIi < liminf^oo H^n^Hi from the weak lower semi- 
continuity of the H 1 -norm (see e.g. |39j [Remark 1.3.1 (iii)] ) , so (I2.27P follows from this, 
(I2.32p . and the asymptotic behavior of (f n ) n and (ip n ,o)n- □ 

We now prove the main result of this section: 

Proposition 2.6. Let T, n and Q, be as in §TT\ and let a G C 2 ([0, T}; H div0 (fi; R). 
Assume that f G W 2 (0, T; V. , H-i), *Po G H 2 , f(0)-H(0)ip G H 2 and f'(0)-H'(0)i/j + 
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i(i7(0)/(0) — H(0) 2 ipo) G %\. Then there is a constant c > depending only on T , Ao, 
A\, A2 and VL, such that the solution tp to (12.201) satisfies: 

lO^/^'Xt)^ < cJ2 [||A fe Vo||i + P'f/^Ww^Ho^)] , J = 0, 1,2, * G [0,71. 

(2.33) 



fc=0 



Proof. The case j = following directly from Lemma I2T51 we first prove (I2.33P for j = 1. 
Taking into account that ^0 G % and / G VT(0,T;'Ho), we have tp G C°([0, 7 1 ]; "H 2 ) H 
C 1 ([0, T]; Ho ) from Lemma [2.21 and ip' is solution to the system 

-i< + i7(t)V = fx ■= f - H'(t)1> in L 2 (0, T; ?{_ 2 ) 

^(o)=i(/(o)- J ff(o)^ b ) I [Z - 6V 

by Lemma E31 Further, / a G W(0, T; K , ft-i) since f G W(0, T; H , H-i), and we 
have 

!|A||l2(0,T;«o) < ||/'||l9(0,T;«o) + ^i7 1 ||^I|l-(0,T;W 1 ) < C 3 (||Vo||l + \\f'\\l?{0,T;Ho)), ( 2 - 35 ) 

according to (12.121) and (12.331) with j = 0, for some constant C3 > depending only on 
T, ^4 0) and Q. This combined with (12.261) entails 

W(t)\\ < c 4 (||Vo||i + II Apollo), t G [0,T], (2.36) 

where C4 > depends on T, Ao, Ai and f2 as well. Moreover, bearing in mind that 
f[{t) = f"{t) - H'(t)ip'(t) - H"(t)ip(t), we find that 

||/i||l2(0,T;«_i) < ||/"||l 2 (0,T;W-i) + hT\\^ ||l°°(0,T;« ) + ^ || tjj || L°°(0,T;Wo) ) 

which, together with flZgfl and fl23oD - fl23oD . proves (I2T331 for j = 1. 

Further, in light of fl2~3l) . A^ e %i and /1 G W(0, T; H , U-i), we know from 
Lemma O that tp G C 1 ([0, T]; Hi). This shows that / a G W(0,T;% ) and, A^ being 
taken in "H 2 , Lemma [2.41 guarantees that ip" is solution to the system 

J + #(^" = / 2 := /" - 2H'{tW - in L 2 (0, T; H- 2 ) (9 ^ 

\ <(0) = i(/'(0) - H'(0)ik) - J7(0)/(0) + J7(0) 2 ^o. 1 J 

Finally applying Lemma I2T21 to (12.341) once more we see that ip e C 1 ([0, T]; "H 2 ) H 
C 2 ([0, T] ; "H ) • Therefore f 2 G W(0, T; ft , H-i), and the case j = 2 follows from (12371) . 
by arguing in the same way as in the derivation of (I2.33P with j = 1, from (I2.35p . □ 
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3 Stability inequality 



Let O, ao and M be the same as in §1.31 Then we consider the time depend magnetic 
Hamiltonian H a (t) = H(t), t G [0, T], defined in §2.21 and associated to a magnetic 
potential of the form a(t,x) := x^a^); where 

a G A(a , M) and x e C 3 ([0, T]; R) is such that x(0) = and x'(0) ^ 0. (3.1) 

Notice in this case that the operators Bj^ a (t) = Bj(t), j = 1,2,3, defined in §2.3[ have 
the following expression: Bj^(t) = x^\t)a..(iV + x(£)a), where x^if) '■= ^r(0- 

3.1 The inverse problem 

In this section we discuss the inverse problem of determining the magnetic potential 
a G A(a , M) from the measurement on a part T + of the boundary T, of Neumann data 
of the solution u (more precisely the first and second order time derivatives of the flux 
of u) to the problem (11.11) . where u is a suitable real- valued function on Q, satisfying: 



In light of (11 .11) and Lemma |2~2| and since Uq G "H 2 , we have u G C°([0, T]; "H 2 ) H 
C 1 ([0, T\ \ %o), and u' is solution to the system 



according to Lemma [2^1 Further, using that fi G W(0, T; H , "H_i) and Au G "Hi, 
it follows from (J3l| and Lemma O that v! G C^T;?^) R C°([0, T]; 7f_i), which, 
in turn, proves that /1 G 1^(0, T; Tin). Thus, taking into account that Au G H 2 , we 
obtain that u G C 1 ([0, T];H 2 ) H C 2 ([0, T]; % ) from (E3D and LemmaO Moreover, by 
Lemma 12.44 M " is solution to the problem 



with / a := -H'l(t)u - 2H'Jt)u' . Finally, by substituting it" (resp. / 2 , A 2 w and (ED) 
for u' (resp. /1, Aw and (13. 3ft ) in the above reasoning, we end up getting that u G 

C 2 ([0,T];^ 2 )nC 3 ([0,T];^ ). 

Further, a and a being taken in A(a , M), let u and u be the respective C 2 ([0, T]; % 2 )n 
C 3 ([0, T}; "H )-solutions to (II. ip and the system 



A k u G U 2 , k = 0,1,2. 



(3.2) 




(3.3) 




(3.4) 




(3.5) 



15 



Then, bearing in mind that V.a = V.a = 0, we find using basic computations that 
v := u — u G Hi satisfies 



-iv'(t, x) + Hjt)v{t, x) = f(t, x), (t, x) G Q% 
v(t,x) = 0, (t,x)e^ 
v(0,x) = 0, x eft, 



(3.6) 



where 



/ := X (a - a). (2iV + *(a + a)) u G W 2 (0, T; ft , ft_x). 



(3.7) 



Since / G W(0,T;ft ), we have v G C°([0, T]; H 2 ) n (^([O, T]; % ) by Lemma Q and 



according to Lemma [2T4l As g G W(0, T; 7^ 0) H-i), (13.81) and Lemma [2721 yield «; G 
C°([0, T}; Hi), which, in turn, implies that g actually belongs to W(0, T; Hq). Therefore, 
we get that w G C°([0, T); H 2 ) H (^([O, T]; Hq) and y := w' is solution to 

f x) + H a (t)y(t, x) = q(t, x) := g'{t, x) - H' a (t)w(t, x), (t, x) G Q+ 



by arguing as before. Finally, taking into account that (a — a).Vuo G Hi (since a and 
a coincide in a neighborhood 7 of Y) and q G W(0, T; % , H-i), we find out from (13. 9p 
and Lemma ETJ1 that y is uniquely defined in C°([0, T]; Hi) n (^([O, T]; %_i). 

Though, in this framework, additional technical difficulties arise from the presence 
of first order spatial differential terms in the expression of H a (t), the method used 
in §3.31 to establish a stability estimate on a is inspired from [10]. This method of 
symmetrization already used in [15] and [H] for Schrodinger systems allows us to center 
the problem around the initial condition which is a given data in our problem and 
avoid the use of another data at time 9 > 0. This method preliminarily requires that 
a Carleman inequality for the solution y to (13. 9ft be established in Qt '■— (— T, T) x Q. 
To this purpose we first extend x m an °dd function on [— T, T], and introduce the 
corresponding operators H a (t) and B a (t) for every t G [— T, T]. Then we extend u (resp. 
u) on := (— T, 0) x Q by setting u(t, x) = u(—t,x) (resp. u(t, x) = u(—t,x)) for 
(t,x) G i n such a way that u (resp. u) satisfies (11.11) (resp. (13. 5p ) in Q T . Since u is 
real-valued by assumption, the mappings t t— > u(t,x) and i i— > u(t,x) are actually both 
continuous at t — for a.e. x G f2. Further, by extending the above definitions of v and 
/ in Qt, we get that v(t,x) = v(—t,x) and f(t,x) = f(—t,x) for all (t,x) G Q^. Hence 
f is solution to (I3.6P in Qt, and, due to the initial condition v(0, x) — 0, t h-> t>(t, x) is 
actually continuous at £ = for a.e. x G f2. Similarly we may define w and g on 



w := v' is thus solution to 




(3.8) 




(t,x) G £+ 



(3.9) 
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from the extended definitions of v and /, as w(t, x) = —w(—t, x) and g(t, x) = —g(—t, x) 
for (t, x) G Qt- This combined with the initial condition w(0,x) = guarantees that 
t i — y w(t,x) is continuous at t = for a.e. x G Q, and it is easy to check that w is 
solution to (13. 8p in Q T . Finally, by setting q(t,x) := q(—t,x) and y(t,x) := y(—t,x) for 
(t, x) G Qt, we end up getting that y is solution to ( 13.91) in Qt- Moreover, in light of the 
identity 2/(0, x) = — 2x'(0)(a(x) — a(x)).Vwo(a;) and the condition u (x) G K imposed 
on u , the mapping t \-> y(t, x) turns out to be continuous at t — 0. 



3.2 Carleman estimate 

In this section we prove a Carleman inequality for the solution y to the system ( 13. 91) . 
extended to Qt- This is a powerful tool in the area of inverse problems which was 
introduced by A. L. Bugkheim and M. V. Klibanov in [20] (see also e.g. [12] and 
[43]). So far, the method defined in [20] is the only one enabling to prove uniqueness 
and stability results for inverse problems with finite measurement data in the multi- 
dimensional spatial case (i.e. n > 2). 

Let us first recall some useful result borrowed from [T3] , establishing a global Carleman 
estimate for any (sufficiently smooth) function q, defined in Qt, which vanishes on 
S T := (— T, T) x dQ, and the Schrodinger operator L acting in (C^)'(Qt), as, 

L:=id t + A. (3.10) 

This preliminarily requires that we introduce an open subset r + of T, together with 
some function (3 G C 4 (fi;M + ), satisfying the following conditions: 

Assumption 3.1. (a) 3C*o > such that we have \V(3(x)\ > Co for all x G fi; 

(b) dj(a) := VP(a).v(a) < for all aeT' := T\T + ; 

(c) 3Ai > 0, 3e > such that we have A|V/3(x).C| 2 + £> 2 /3(C,C) > e|C| 2 for all ( G 



and X > Ai, where D 2 f3 = ( «f 2 £, ) , D 2 (3((,() denotes the C n -scalar 

product of D 2 j3( with (. 



Notice from Assumption 13.1( c) that /3 is pseudo-convex with respect to the operator 

—A. We refer to [15] [Formula (6)] for actual examples of an open subset r + and a 
weight /3, fulfilling Assumption 13.11 

Further we put 

j3 := (3 + K, where K := mH/^Hoo for some m > 1, (3-11) 
and then define for every A > the two following weight functions: 

e X/3(x) e 2XK _ e X/3(x) 

(p(t,x) = — — r and ri(t, x) = — r-— r, (t,x) G Qt- (3.12) 

^ v ' (T + t)(T-t) /v ' (T + t)(T-ty v ' y ' 
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Finally, following the idea of [15] [Proposition 3], we introduce two operators acting in 
(CgWr), 

Mi := id t + A + s 2 |Vn| 2 and M 2 := isr/ + 2sVn.V + s(An), (3.13) 

in such a way that M\ + M 2 = e~ sv Le sv , where L is given by (13. 10p . Then, by arguing 
in the exact same way as in the derivation of [15] [Proposition 1], we obtain the: 

Proposition 3.2. Let T , n and Q be as in let (3 be given by (13.111) . where T + C T 
and (3 G C 4 (f2;IR + ) fulfill Assumption \3.1\ let tp and rj be as in (I3.12p . and let L, M\ 
and M 2 be defined by H3. 1Q[) - A3. 13[) . Then there exist three constants Xq > 0, s > 0, 
and Co = Cq(T, Q, T + , X , sq) > 0, such that we have 

I(q) < C ^sxj^l^e- 2s ^^(t,a)d u ^a)\dMt,^ 

'(3.14) 

for all X > Aoj all s > So, and all q G L 2 (— T, T'Tii) satisfying Lq G L 2 (Qt) and 
d v q G L 2 (-T,T;L 2 (r)), where 

I(q) := S 3 A 4 ||e-V /2 9||^( QT ) + S A||e-V /2 |V g |||^ (QT) +^ WM^qWh^y (3-15) 

j=l,2 



We turn now to proving a Carleman estimate for the solution y to the system (13. 9 j) 
in Q T , i.e. 

-iy'(t, x) + H a (t)y(t, x) = q(t, x), (t, x) G Q T 

y(t, x) = 0, (t, x) G S T := (-T, T) x T (3.16) 

y(0,x) = -2x'(0)(a(x) - a(x)).Vu (x), x G O, 

with the aid of Proposition 13.21 

Proposition 3.3. Let n and Q be as in Lemma \2. H /ei /3, tp and n be the same as in 
Proposition Iff.ffi let x fulfill (13. ip . let uq satisfy (13 ,2ft . Zet a and a belong to A(ao, M), 
and /ei /(y) fre defined by ( 13. 151) . where y denotes the solution to (I3.16p . Taen there 
exist three constant X > 0, s > 0, and C\ = Ci(T, M, Q, T + , A , s , u ) > 0, such that 
we have 

I(y) < Ci (sX f f e- 2s ^v(t,a)d u {3(a)\d u p(t,a)\ 2 dadt 

^5 _ r^ll 2 . 



+ e~ s " a- a 



lL 2 (Q T ) n y ' 



for any X > Aq and any s > s . 
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Proof. In light of (I3.10P we see that the first equation in (13.81) (extended to Qt) reads 
Lw = 2%a.(iV + x a ) w ~ 9-, the r.h.s. of this equality being in L 2 (Qt)- From this and 
Proposition 13.21 then follows that 

I(w) < C (sX J J e- 2sr >^<p(t,(r)\d v w(t,a)\ 2 d v p(t,a) dadt 

+ \\e~ sri f\\h { Q T ) + E + I V p| 2 ) 1/2 |Il 2( q t) ) , (3-17) 

p=V,W / 

provided A and s are taken sufficiently large. Similarly, we deduce from ( 13 . 1 6[) that 
I(y) < C (sX j^e- 2s ^^a)\d u y{t,a)\ 2 d u (3(t,a) dadt 



L 2 (Qt) " r 

P 



E H^dPp + lVpl 2 ) 1 / 2 !!^), (3.18) 

=v,w,y / 



the constant Co being the same as in (I3.17p . 

Further, by applying Lemma [374| whose proof is postponed to the end of §3.2[ succes- 
sively to v(x, t) = J '«;((,i)d( and Vv(x,t) = J" * Vu>(£, x)d£, we find out for all A > A 
and all s > 0, that 

He- (H 2 + \Vv\ 2 y/ 2 \\ 2 LHQT) < l\\e-* (\w\ 2 + \Vw\ 2 fl 2 \\l HQTy 

where the constant k = k(T, Xq) > depends only on T and Ao- Therefore, upon 
choosing s sufficiently large and eventually substituting 2C for C , we may actually 
remove v from the sum in the r.h.s. of both (13.171) and (I3.18p . Moreover ||e _,sr? (|p| 2 + 
| Vp| 2 ) 1 / 2 || 2 2( -Q T - ) , for p = w,y, being made arbitrarily small w.r.t. I(p), by taking s 
large enough, according to (13. lip - ( 137121) and ( 13. 151) . we may rewrite (13.171) and ( 13. 18ft 
as, respectively, 

I(w) < C 2 (sxJ T T J^e- 2s ^cp(t,a)\d u w(t^)\ 2 d u (3(t,a) dadt + \\e~ sri f'\\l HQT ^ , 
and 

I(y) < C 2 {^sX J T ^J^e- 2s ^^a)\d v y{t,a)\ 2 d v P{t,a)dadt 

+\\e- s y"\\l HQT) + i(w] 

Finally the result follows from this, the two following expressions, 

/' = (a - a).[ X (a + a)(2 x '« + X u) + 2i{ X 'Vu + X Vu% 
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and 



/" = (a - a).[2i( X "W + 2 X 'Vu' + X Vu") + (a + a)(2 X ' 2 + XX ")u + A XX 'u' + X 2 u' '}, 

arising from (13.71) by direct calculation, and Proposition ^. 61 (in the particular case where 
the source term / = 0). □ 

To complete the proof of Proposition 13.3} it remains to establish Lemma 13.41 Its 
proof, although very similar to the one of jH] [Lemma 2.1] or |45j [Lemma 3.1.1], is 
detailed below for the convenience of the reader. 

Lemma 3.4. Let T , n and Q be as in and let r\ be defined by (13.121) . Then for all 
Ao > ; there exists a constant k = k(T, Ao) > depending only on T and \q, such that 
we have 

„T „ „t 2 

K , 



-2sr/(t,x) 



T Jn 



dxdt < -||e >||l2( Qt) , 



for every p G L 2 (Q T ), all A > A and all s > 0. 
Proof. In light of (l3TTT) -( l3~T2l) . it holds true that 



d t r)(t,x) 



2ta{x) 



(T + t) 2 (T-ty 



with a(x) = e 2XK - e A/3(x) > a := e 2X ° K - e XaK/m > 



(3.19) 

for every (t : x) € Qt and A > A . Further, for all 5 G (0, T), the integral Is :— 

2 

J-tIs fa e~ 2sr,{t ' x) fo p{£, x)d£ dxdt satisfying 



Is < 



[ T+S [ e~ 2s ^ x h ( f \p(Z,x)\ 2 dt) dxdt, 
J-t+s Jn \Jo J 



from the Cauchy-Schwarz inequality, we get from (13 . 1 9[) that 




T-S 



e- 2sr ' {t ' x) d t r](t, x) ( / \p{t x)\ 2 di ) dtdx. 



2T 4 

"o Jn J-T+6 

Moreover, for a.e x G Q, an integration by parts show that 



(3.20) 



T-S 



T+S 



e- 2s ^d t r)(t,x)( / \p(t,x)\ 2 d{ )dt 



1 

2s 



T-S 



T+S 



-2sr/(t,x) 



\p(t, x)\ 2 dt - p(T -S,x)+ p(-T + 5, x) 



where p(t,x) : = e" 2 ^*'^ ( /* \p(£,x)\ 2 d£). This together with fl3T2TJ|) yields 



Is < 



rp4 



T-S 



-2sr](t,x) 



\p{t,x)\ 2 dxdt 



«oS \J-t+s Jn 

-2(T - 5) I (p(T -5,x)+ p(-T + 5, x))dx 



(3.21) 



20 



Since f n p(±(T - S),x)dx < e" 2sao/(5(2T - 5)) |H|^_ TjT;C) by fl3l2|) and ( l3~T9j) . we get 
that 

lim(T -5) [ (p(T -6,x)+ p(-T + 5, x))dx = 0, 
54.0 y n 

and the result follows from this, by taking the limit as 8 J, in f)3.2ip . □ 

3.3 Global Lipschitz stability inequality: proof of Theorem HJj 

In this section we establish the Lipschitz-type stability inequality for admissible magnetic 
potential vectors a £ A(a ,M), stated in Theorem 11.11 The proof essentially relies on 
the following: 

Lemma 3.5. Under the assumptions of Propositio ns. 31 lety denote the C°([— T, T]; Tii) - 
solution to (13.161) and put X := \\e~ sv ^°'^y(0, -)IIl 2 (q)- Then there are three constants 
A > 0, so > 0, and C3 = C(T, M, Q, T + , X , sq, uo) > 0, such that we have 

+s- 1 A- 1 ||e- sr ' (0 " ) (a-a) 



lL 2 (Q) 

for all A > A and s > so- 

Proof. Set if) := e~ sri y. Bearing in mind that ip(—T, .) = 0, we find out that 



X = 
whence 

X = 21m 



Or / rO 



TJfl \J-TJQ 



d t \4>{t,x)\ 2 dxdt = 2Re ( / / $'(t, x)ip(t, x)dxdt) , 



t Ja 




(iip'{t, x) + Aip(t, x) + s 2 \Vr](t, x)\ 2 4){t, x)) ip{t, x)dxdt J 

2Im^y J Mirp(t,x)^(t,x)dtdx S j . 

Therefore \X\ < 2||Mi'0||l2(q t )||'0||l2(q t ) from the Cauchy-Schwarz inequality. Since 
inf {(p(t, x), (t,x) £ Qt} > by ( 13 .121) . there is thus a constant C4 > depending only 
on T and Ao, such that we have 

\X\ < C,s-V 2 \- 2 (s 3 \*\\e-^ 2 y\\l {QT) + WM^y)^^) , 

for all s > sq and A > Aq. This, combined with ( 13 . 1 5f) and Proposition 13.31 yields 



1x1 < c,,-v^( E /;/ 

\p=w,y 



e- 2s ^<p(t, a)d v p{t, a)\d v p(t, a) \ 2 dadt 

r+ 



+s- 1 A- 1 ||e- s "(a-a)||2 2(QTr ,, 
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the constant C5 > depending on Q, T + , T, M, X , so and uq. Now the result follows 
immediately from this since rj(t,x) > rj(0,x) for all (t,x) G Qt by (I3.12p . □ 

Armed with Lemma l3"3l we turn now to proving Theorem ll.il Since X = 4x'(0) 2 ||e _ ' sr '( '^(a— 
a).Vn ||o from (13. 9p . Lemma [331 assures us that there are three constants Ao > 0, s > 
and Cq > 0, such that we have 



G 6 (||e-^(a-a).V«o||^ (n) -^A- 2 ||e^°'0(a-a)||^ (n ^ 
< s~ ll2 \- 1 V / / e- 2sv{t ^ip(t,a)d u (3(t,a)\d u p(t,a)\ 2 dadt, (3.22) 



p=w,y 



for all s > Sq and A > Ao- In light of the definitions of w, y, and those of Uqj, Uj, Uj, 
for j = 1, . . . , n, (I3.22p entails 

C 6 (||e-^°'-)(a^^ 
< s-^X' 1 J J e- 2sr >^tp(t,a)d v p(t,a) (j^ l<^(% - «j)(^)| 2 j dadt. 

Summing up the above estimate over j — 1, . . . , n, we get that 
C 6 (\\e~ sr >^DUo(a - a)||^ (n) „ - ns^A-le""^^ - a) 



|2 

lL 2 (f2)' 



j — 1, . . . , n 

h = 1, 2 



dadt, 



(3.23) 

where the n x n real matrix DUq(x), for x G f2, is the same as in (11.21) . Notice that we 
have 

\\DU (x)v\\ C n > ^i(x)\\v\\ C n, iGfl, v e C\ (3.24) 

where { / u :; (x)}™ =1 C i?" denotes the non decreasing sequence of the singular values 
of DUq(x), and ||v||c" stands for the Euclidian norm of v. Moreover each Uqj, j = 
l,...,n, being taken in H p (f2) with p > n/2 + 1, it holds true that uo,j G C 1 (f2), 
whence \i\ G C°(f2;IR + ) from [4"6] [Theorem II. 5.1]. This, combined with (II. 2p . yields 
H := inf^gn /ii(x) > 0, the constant fi depending on Q and {m j}" =1 only. As a conse- 
quence we have ||e- s ^ '-)DC/o(a-a)|| L 2(o)n > ^||e- s ^ ")(a-a)|| L2(n )n by (ET23)1 . and The- 
orem [TTTJ follows directly from this and (I3.23P by choosing s so large that us~ 3 ^ 2 Xq 2 < /i 2 . 



22 



Acknowledgement (s) 



The authors are thankful to M. Bellassoued for fruitful discussions concerning this prob- 
lem and valuable comments. 



References 

G. ESKIN, A new approach to hyperbolic inverse problems, Inverse Problems 22 
(2006), 815-833. 

G. Eskin, Inverse boundary value problems and the Aharonov-Bohm effect, In- 
verse Problems 19 (2003), 49-62. 

Z. Sun, An inverse boundary problem boundary value problem for the Schrddinger 
operator with vector potentials, Trans. Am. Math. Soc. 338 (1992), 953-969. 

G. Nakamura, Z. Sun, G. Uhlmann, Global identif ability for an inverse prob- 
lem , Math. Ann. 303 (1995), 377-388. 

C. F. Tolmasky, Exponentially growing solutions for nonsmooth first-order per- 
turbations of the Laplacian, SIAM J. Math. Anal. 29 no. 1 (1998), 116-133. 

M. Salo, Inverse problem for nonsmooth first order perturbation of the Laplacian, 
Ann. Acad. Sci. Fenn. Math. Dissertations 139 (2004). 

A. PANCHENKO, An inverse problem for the magnetic Schrddinger equation and 
quasi- exponential solutions of nonsmooth partial differential equations, Inverse 
Problems 18 (2002), 1421-1434. 

M. Bellassoued, M. Choulli, Stability estimate for an inverse problem for 
the magnetic Schrddinger equation from the Dirichlet-to-Neumann map, Journal 
of Functional Analysis, 91 no. 258 (2010), 161-195. 

G. Eskin, Inverse problems for the Schrddinger operators with electromagnetic 
potentials in domains with obstacles, Inverse Problems 19 (2003), 985-996. 

M. Salo, Semiclassical pseudodifferential calculus and the reconstruction of a 
magnetic field, Comm. Part. Diff. Equ. 31 no. 11 (2006), 1639-1666. 

G. Eskin, Inverse problems for the Schrddinger equations with time- dependent 
electromagnetic potentials and the Aharonov-Bohm effect, J. Math. Phys. 49 no. 
2 (2008), 022105, 18 pp. 

D. Dos Santos Ferreira, C. E. Kenig, J. Sjostrand, G. Uhlmann, 
Determining a magnetic Schrddinger operator from partial Cauchy data, Comm. 
Math. Phys. 2 (2007), 467-488. 

L. TZOU, Stability estimates for coefficients of magnetic Schrddinger equation 
from full and partial boundary measurements, Comm. Part. Diff. Equ. 33 (2008), 
1911-1952. 



23 



H. Ben JOUD, A stability estimate for an inverse problem for the Schrddinger 
equation in a magnetic field from partial boundary measurements, Inverse Prob- 
lems 25 no. 4 (2009), 45012-45034. 

L. Baudouin, J. -P. Puel, Uniqueness and stability in an inverse problem for 
the Schrddinger equation, Inverse Problems 18 (2002), 1537-1554. 

A. Mercado, A. Osses, L. Rosier, Inverse problems for the Schrddinger equa- 
tion via Carleman inequalities with degenerate weights, Inverse Problems 24 no. 
1 (2008), 015017. 

M. Bellassoued, M. Choulli, Logarithmic stability in the dynamical inverse 
problem for the Schrddinger equation by arbitrary boudary observation, J. Math. 
Pures Appl. 91 no. 3 (2009), 233-255. 

P. Albano, Carleman estimates for the Euler- Bernoulli plate operator, Elec. J. 
of Diff. Equ. 53 (2000), 1-13. 

D. Tataru, Carleman Estimates, Unique Continuation and Controllability for 
anisotropic PDEs, Cont. Math. 209 (1997), 267-279. 

A.L. Bukhgeim, M.V. Klibanov, Uniqueness in the large of a class of multi- 
dimensional inverse problems, Soviet Math. Dokl. 17 (1981), 244-247. 

O. Yu. IMANUVILOV, M. YAMAMOTO, Lipshitz stability in inverse parabolic 
problems by Carleman estimate, Inverse Problems 14 (1998), 1229-1249. 

M. Choulli, M. Yamamoto, Generic well-posedness of a linear inverse 
parabolic problem with respect to diffusion parameters J. Inv. Ill-Posed Problems 
7 no. 3 (1999), 241-254. 

V. Isakov, M. Yamamoto, Carleman estimates with the Neumann boundary 
condition and its applications to the observability inequality and inverse problems, 
Cont. Math. 268 (2000), 191-225. 

A. L. Bukhgeim, J. Cheng, V. Isakov , M. Yamamoto, Uniqueness in 
determining damping coefficients in hyperbolic equations, Analytic Extension For- 
mulas and their Applications (2001), 27-46. 

O. Yu. Imanuvilov, On Carleman estimates for hyperbolic equations Asymp- 
totic Anal. 32 (2002), 185-220. 

M. Bellassoued, Uniqueness and stability in determining the speed of propa- 
gation of second-order hyperbolic equation with variable coefficients, Applicable 
Analysis 83 (2004), 983-1014. 

M. Choulli, M. Yamamoto, Conditional stability in determining a heat source, 
J. Inv. Ill-Posed Problems 12 no. 3 (2004), 233-243. 

M. Bellassoued and M. Yamamoto, Logarithmic stability in determination 
of a coefficient in an acoustic equation by arbitrary boundary observation , J. Math. 
Pures Appl. 85 (2006), 193-224. 



24 



M. Cristofol, P. Gaitan, H. Ramoul, Inverse problems for a 2X2 reaction- 
diffusion system using a Carleman estimate with one observation, Inverse Prob- 
lems, 22 , (2006), 1561-1573. 

M. V. KLIBANOV, M. YAMAMOTO, Lipschitz stability of an inverse problem for 
an accoustic equation, Applicable Analysis 85 (2006), 515-538. 

Bellassoued and M. Yamamoto, Lipschitz stability in determining density 
and two Lame's coefficients , J. Math. Anal. Appl. 329 (2007), 1240-1259. 

M. Cristofol, L. Roques, Biological invasions: deriving the region at risk from 
partial measurements, Math. Biosci. 215 no. 2 (2008), 158-166. 

A. Benabdallah, M. Cristofol, P. Gaitan, M. Yamamoto, Inverse prob- 
lem for a parabolic system with two components by measurements of one compo- 
nent, Applicable Analysis 88 no. 5 (2009), 683-710. 

H. Ammari, G. Uhlmann, Reconstruction of the potential from partial Cauchy 
data for the Schrddinger equation, Indiana Univ. Math. J. 53 (2004), 169-183. 

B. Helffer, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, M. P. 
Owen, Nodal sets for the groundstate of the Schrddinger operator with zero mag- 
netic field in a non simply connected domain, Comm. Math. Phys. 202 (1999), 
629-649. 

R. Dautray, J.-L. Lions, Mathematical analysis and Numerical Methods for 
Science and Technology, Vol. 5, Evolution problems, Springer Verlag, Berlin Hei- 
delberg, 1992. 

J. Li, M. Yamamoto, J. Zou, Conditional stability and numerical reconstruc- 
tion of initial temperature, Comm. Pure. Appl. Anal 8 (2009), 361-382. 

J.-L. Lions, E. Magenes, Problemes aux limites non homogenes et applications, 
vol. 1 & 2, Dunod (1968). 

T. Cazenave, Semi-linear Schrddinger equations, Courant Lecture Notes, Vol. 
10, AMS & Courant Institute of Mathematical Sciences, 2003. 

O. Yu. Imanuvilov, M. Yamamoto, Global Lipshitz stability in an inverse 
hyperbolic problem by interior observations, Inverse Problems 17 no. 4 (2001), 
717-728. 

L. Cardoulis, M. Cristofol, P. Gaitan, Inverse problems for the 
Schrddinger operator in an unbounded strip, J. Inv. Ill-Posed Problems 16 no. 
2 (2008), 127-146. 

A.L. Bukhgeim, Volterra Equations and Inverse Problems, Inverse and Ill-Posed 
Problems Series, VSP, Utrecht (1999). 

M. V. Klibanov, Inverse problems and Carleman estimates, Inverse Problems 
8 (1992), 575-596. 



25 



[44] M. V. Klibanov, Global uniqueness of a multidimensionnal inverse problem for 
a non linear parabolic equation by a Carleman estimate, Inverse Problems 20 
(2004), 1003-1032. 

[45] M. V. Klibanov, A. A. Timonov, Carleman Estimates for Coefficient Inverse 
Problems and Numerical Applications, VSP, Utrecht (2004). 

[46] T. Kato, Perturbation Theory for Linear Operators, Die Grundlehren der math- 
ematischen Wissenschaften 132, Springer- Verlag New York, Inc., New York 1966. 



26 



